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On Evaluating {*2e"°"°” dx by Contour 
Integration Round a Parallelogram 


Darrell Desbrow 


Introduction. The usefulness of contour integration for the evaluation of infinite 
integrals can hardly be over-stressed. Strangely, however, and this can be some- 
thing of an anti-climax for the student, the method finds no ready application to 
what is arguably the most useful and commonplace such integral, viz., the probabil- 
ity integral 


+ 00 

J ede = Vr. (1) 
Indeed it was widely held at one time that contour integration could not achieve 
this evaluation. By now, however, it is fairly well known that the probability 
integral can indeed be evaluated by a contour integration of a suitable integrand 
around a suitable skew (i.e., non-rectangular) parallelogram. Such a proof was 
given by Mirsky [3] as long ago as 1949. There is an account of such an evaluation 
in [2]. It will be seen from either reference that neither the integrand nor the 
parallelogram contour employed is at all evident. 

The related Fresnel integrals 


2 us stato 
J cos x° dx = =| sin x* dx, 
ae p) = 


which can be written in complex form as 


foe dx = i =71)a (2) 


are easier to handle and can be evaluated by integration of a suitable integrand 
round a suitable standard rectangle [1]. Again, however, the integrand in particular 
is not at all evident. 

For the student familiar with the notion that integration around rectangles, and 
so by extension around parallelograms, is a useful technique, the difficulty still 
remains, regarding the probability and Fresnel integrals in particular, of finding 
suitable integrands and contours in order to make the desired evaluations. It 
hardly satisfies, even if it suffices, to be bidden without explanation to integrate 
this function round that parallelogram. And even if one is, the doubt lingers 
whether it is really necessary to have to resort to a skew parallelogram for the 
probability integral when a rectangle serves so well for the Fresnel integrals. The 
purpose here is to shed light on these matters but in relation to the more general 
integral 


+ co 
J ere 2 Nye: (3) 


—e* 


with a,b € C suitably restricted for convergence, of which the probability and 
Fresnel integrals are the respective special cases b := 0,a := —1, —i. We analyse 
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the two cases 


i) rea < 0, into which the probability integral (1) falls, and 
ii) rea = 0 > ima and im b < 0, into which the Fresnel integral (2) falls. 


It transpires that, in these cases, the integral (3) evaluates to 

+i ied ea 

ar 
wherein va is the principal square root of a and the sign is fixed by a definite rule. 
In itself this evaluation is a secondary issue; being a rather routine, if technical, 
matter once an integrand and contour have been specified. There are, moreover, 
other evaluations for the probability and Fresnel integrals easier than the contour 
integrations proposed here. The primary aim here is to demonstrate to the student, 
as regards the probability and Fresnel integrals in particular, that an evaluation by 
contour integration is possible, how an integrand and parallelogram contour 
suitable for the evaluation can be hit upon and to show why, given the tack taken, 
rectangular contours cannot work for the probability integral in particular. 

A third case, viz., iii) rea = 0 < ima and im b > 0, can be reduced to the case 

ii) by conjugation, on noting that 


pe (x-2b) — a@x(x—2b) 
i ere dx = [ e** dx, 


—@ 


/ 
when either integral exists. By this ruse of conjugation we may, and hereafter do, 
assume that im a < 0 when rea = 0. 

Throughout, for any z © C, yz denotes the (principal) square root of z, viz., 
vre®/? when z =: re’’ with r> 0, —a < 60 <7. The square roots of z are then 
precisely + yz. 


Determination of integrand and parallelogram. Let a,b <€ C with a # 0. For any 
R, with 0 < R < +, let 


I(R) a= ie eax(x—2b) dx. 


In case R # +~%, let y be the parallelogram ABCD with vertices +a + RB as 
shown in Figure 1, where a € R, with a > 0 and B € C, with 0 < arg B< a. 


Figure 1. Parallelogram contour y. 
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We consider the contour integral J := |, f(z) dz, wherein f and y are to be 
chosen so as to permit evaluation of /7; bel er 20) de = I(+), as the Cauchy 
principal value lim, _, ,../(R). 

Towards the evaluation of J, integration along the sides AB and CD of y, 
contributes 


Bf [f(#B + @) —f(xB - @)] a 


which is JCR) for the choice B[f(xB + a) — f(xB — a)] = exp ax(x — 2b). On 
writing z=xB+a, 


fat 9n) = sex a — =\(= c - 26] re 


B B 
so that 
ne et 2b 4 
(Z)= poe 2? a)(zZ—a B) (4) 
= h(0)exp pele Sa bpy, (5) 


This suggests that a suitable integrand might be got by solving the equation 
f(z) — f(z — 2a) = A(z) for f(z) and to that end we have the following proposi- 
tion, whose proof is by direct substitution. 


Proposition 1. Suppose that h(z — 2a) = h(z)g(z) for some function g with period 
2a and each z € UC. Then 


A) 
f(z )= 7 1 1—g(z) 
solves f(z) — f(z — 2a) = h(z), when g(z) ¥ 1. 


We now seek to determine a suitable function g. From (4) and (5), 
a 
h(z — 2a) =h(0)exp Bi? — 2a)(z — 4a — 2bB) 


Z 


= h(z)exp —5-e 


4aa 
B xp “pee 


Thus 


8a’a 4aa 
Be exp "Bz ee — z) 


serves as a suitable g provided it has period 2a, thus provided 2a = 27iB*/4aa 
or 4a*a = wif’. Since this equation merely determines the ratio B/a, we may (at 
the expense of replacing B by B/2qa) assume that a:= 1/2. With that choice, 
a = wi” so that 


exp 


va 
p= ie, (6) 
where w := vi = e'7/* and 
the sign is chosen to ensure that 0 < arg B < 7w. (7) 
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Note from the equation a = 7wif* that B is real if and only if rea = 0 < ima. 
Thus £ is not real, by assumption already made that ima < 0 when rea = 0. 
Then with these choices for a and B, 


1 
h(z) = aoe mi(z — 1/2)(z — 1/2 — 2bB) (8) 


= h(O)expziz(z — 1 — 2bB), 
g(z) = exp2mi(bB —z), and 


f(z) = 
2 a 
1 — g(z) 
It remains to verify that these choices yield the required evaluation. 

Before that, however, observe that the contour chosen is rectangular if and only 
if B takes the form iA for some A > 0. From the equation a = wif’, that is the 
case if and only if rea = 0 > im a; a condition satisfied by the Fresnel integral but 
not the probability integral. Thus, given this approach, a skew parallelogram 
contour is inevitable for evaluation of the probability integral. 


Evaluation of the integral. Clearly the singularities of f are simple poles at z = bB 
(mod 1), where g(z) = 1. [We leave aside the case when one of these poles falls on 
the side AB of y, so that a second falls on the side CD. This is the case, for large 
R at least, precisely when 2im b is an odd integral multiple of im(1/8). Then we 
must make small indentations in the usual way or make some other provision in 
the analysis, e.g., by shifting y half a unit to the right.] Outside that exceptional 
case, just one of the poles, say bB — k for some integer k, is inside the parallelo- 
gram when R > |b|. Its residue there is 
h(bB —k) 7 1 ACB k) _ hb ay 


since g(bB —n) =1 for each n € Z. Further, since h(z-— 1) = hA(z)g(z) for 
each z € C, it follows by easy inductions that h(bB — n) = h(bB) for each n € Z. 
Furthermore from (6) and (8), one checks that h(bB) = w e~*"'/B. Hence by the 
residue theorem, 


I=2 Ab ee ey) eae 
eg vse im) ea) ae 


On the other hand, 

Poe = 

AB CD BC DA 
1/2 1/2 

I(R) + [ f(RB-x)de+ [~ f(-RB +x) dx 
—1/2 =172 


:I(R) +J(R) + K(R), | say. 
We shall show that, for suitable restriction upon a and b, J(R) > 0as R > +, 


A similar analysis shows likewise that KCR) — 0 as R ~ +. Thus, proceeding to 
the limit as R — +, we conclude that 


IT — 


f(z) dz+ f(z) ad 


—@W 


+ 00 W Vr 
if e@X(x—-25) dy = lim I(R) =I = ~e ab* ee a ae) (9) 
Ro +0 Bp a 


with the sign chosen as in (7). This is the desired evaluation. 
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It remains to verify that J(R) > 0 as R > +, 


i h(RB- x) 


|J(R)| = ipl epee b OR BSa 


[o f(RB~ x) ds - 
20/2 
2 ~ |h( RB — x)| 


-12 |1—e"| 
where for |x| < 1/2, 
w(x) = 27i[(b + R)B —x| = 27(b + R) Bi — 2axi. (11) 


In order to proceed with the estimate, we use the 


dx, (10) 


Proposition 2. For any w=u+iveEC, |l —e”|>1/2 when |u| is sufficiently 
large. 


1 as u—> —%, 
Now from (11), as R > +, 
|u(x)| =|re w(x)| = 2ajim [(b + R)B]| > +2. 
Note that this occurs independently of x. Thus for all sufficiently large R, 


independently of x, [1 — e”|> 1/2 for |x| < 1/2. Hence from (10), for all 
sufficiently large R, 


J(R)| <2/ 


Seley 


“ |a(RB ~ x) | dr 


- 2n(O)If exp mi(RB —x)(RB—x-—1-2bB)|dx 


1/2 


= 2|A(0)|f “ e% de, (12) 


7? 
where 
(x) = re [wi(RB — x)(RB-x - 1 — 2bB)] 
= (rea)R* + c(x)R+ d(x), (13) 
for some real-valued continuous functions c,d. Thus for |x| < 1/2, 
b(x) < (rea)R*?+CR+D =: 9(R), say, 
where C:= max c(x) and D:= max d(x). 
|x|<1/2 Ix|<1/2 


Case i): rea <0. As R > +, o(R) > —©&. Further from (12), for all sufficiently 
large R, 


|J(R)| < 2|A(0) yp e dy < 2|h(0) ler. 
=e Y 


Thus J(R) > 0as R-> +”, ca 
Case ii): re a=0> ima and imb <0. Then a =: —iA, for some A > 0. By 
choice of 8B, —iA =a = wif’, and arg B > 0, so that 
vA — 
= ji—. 
a (14) 
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Now for |x| < 1/2, from (13) with rea = 0, 6(x) = c(x)R + d(x) < c(x)R + D, 
where 


c(x) 


re [—i(2Bx + B + 2bB”)| = re|-—aiB(2x + 1) — 27B7ibd| 
—VAn(2x+1)+2Aim), from (14). 
From (12), for all sufficiently large R, 


LJ(R)| < 2|A(0)| f°" 0 ax 
ay? 


<|h(0)| e248 im 9) +D ie o- VATROXx+1) dy 
= 1/9 
(1 = eroise®) 2AR im b 


=|h(0) jeraem b)+D = K(1 -_ en 2VAmR) “ . 


2VAT R 


where K = |h(0)|e? /2V Az. Both e~* 47" and e74*'™°/R tend to0 as R > +. 
Thus J(R) > 0as R > 4+, fe 


In this second case note, from (9), that 


[et dy = © ena? a ea 

a B VA 

Evaluation of the probability and Fresnel integrals. Let b := 0. The respective 
particular choices a := —1 and a := —i yield from (6) 


a d te 
oa ae oo ee sae 


on choosing the sign to ensure that arg B > 0 in each case. These yield respec- 
tively, from (9), 


+o 4 
e the probability integral i e* dx= ae = ya, and 


+0 (7 7 
¢ the Fresnel integral [ 6 a= ra = -iwV7 = \ ry (1 —i). 


Note that, in the case of the probability integral, the parallelogram contour is 
skewed at an angle of 77/4 from the vertical. In the cases of the Fresnel integral 
the contour is a rectangle in standard position. 
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